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We present measurements of the branching fraction and CP-violating asymmetries for neutral 
B decays to D* ± D^ 1 . The measurement uses a data sample of approximately 88 million T(4«S) 
— > BB decays collected with the BABAR detector at the SLAC PEP-II asymmetric-energy e + -e~ 
collider. By fully reconstructing the D* ± D^ decay products, we measure the branching fraction 
to be (8.8 ± 1.0 ± 1.3) x 10 -4 and the time-integrated CP-violating asymmetry between the rates 
to D*~D + and D* + D~ to be A = -0.03 ± 0.11 ± 0.05. We also measure the time-dependent CP- 
violating asymmetry parameters to be S-+ = -0.24 ± 0.69 ± 0.12, C_+ = -0.22 ± 0.37 ± 0.10 for 
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B -> D*~D + and S+- = -0.82 ± 0.75 ± 0.14, C+_ = -0.47 ± 0.40 ± 0.12 for B 
case the first error is statistical and the second error is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



D* + D~. In each 



Within the Standard Model (SM) of electroweak inter- 
actions, CP violation is the result of a complex phase 
in V, the Cabbibo-Kobayashi-Maskawa (CKM) quark 
mixing matrix pj. In the SM, the time-dependent CP- 
violating asymmetries in B — > D* ± £) T decays are re- 
lated to the angle j3 = arg[— V^dV^/VtdV^] . We present 
a measurement of the branching fraction and a first mea- 
surement of CP -violating asymmetries in B — ► £)* ± P) =F 
decays using a sample of 87.9 ± 1.0 million BB decays. 

As recent measurements of the parameter sin2/3 using 
the quark process b — > ccs have shown, CP is violated 
in the neutral B- meson system 0, @- The measured 
asymmetries are currently consistent with the SM expec- 
tation 0. In order to search for additional sources of 
CP violation from new physics processes, different quark 
decays such as b — > ccd must be examined. 

In b — > ccd processes (for example, B — > D* ± _D T de- 
cays; see Fig. 1), penguin contributions containing a dif- 
ferent weak phase than the tree are not expected to be 
as highly suppressed as in b — > ccs decays; thus the rela- 
tion of the time-dependent CP- violating asymmetries in 
b — ► ccd decays to /3 is less exact than in decays such as 
B° — > J/ij)K®. However, the contribution from additional 
weak phases in time-dependent asymmetries in b — > ccd 
due to purely SM processes is still expected to be fairly 
small, of order A/3 = 0.1 in a simplified model [1,0. A 
variety of beyond-SM processes, which can provide addi- 
tional sources of CP violation, can greatly increase this 
contribution, up to A/3 w 0.6 in some models 0. 
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FIG. 1: The leading-order Feynman graphs for B° — * D* ± D^ 
decay: a) tree diagram and b) penguin diagram. 

CP- violating asymmetries in B — ► D* ± £) T are due to 
interference between the decay amplitudes and the B°P° 
mixing amplitude, as well as interference between tree 
and penguin decay amplitudes. The decay rate distribu- 
tions f^, where the superscript +(— ) refers to whether 
the tagging meson, Btag, was P° (P°), are given by 

e -|At|/T 

J At 

[1 ± S sm(Am d At) t- Ccos(Am d Ai)], (1) 



where r is the mean B° lifetime, Am^ is the B°B° mix- 
ing frequency, and At = t leco — itag is the time elapsed 
between the B decay to D* ± D =F and the decay of Btag- 
Separate S and C parameters are fitted for the two decays 
D*-D+ and B>*+B>-, resulting in the four fitted CP- 
violation parameters {<S_ + , C_ + , S + _, C + _}. The time- 
integrated asymmetry A between the rates to D*~D + 
and D* + D~ is defined as 



A 



N 



D' + D- 



D—D-i 



N D , +D - +N D ,- D , 



(2) 



The states D*~ D+ and D*+ D~ are not CP eigen- 
states. The formalism of time evolution for non-CP 
eigenstate vector-pseudoscalar decays is given in Ref . . 
In the case of equal amplitudes for B — * D*~D + and 
B — > D* + D~, one expects that at tree level C_ + = 
C + _ = and S_ + = S + _ = - sin2/3. 

A detailed description of the BABAR detector is pre- 
sented in Ref. 0- Charged particle momenta are 
measured in a tracking system consisting of a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
hexagonal-cell wire drift chamber (DCH) filled with a gas 
mixture of helium and isobutane. The SVT and DCH 
operate within a 1.5 T solenoidal field. Photons are de- 
tected and their energies are measured in a CsI(Tl) elec- 
tromagnetic calorimeter. Muons are identified in the in- 
strumented flux return (IFR), composed of resistive plate 
chambers and layers of iron that return the magnetic 
flux of the solenoid. A detector of internally reflected 
Cherenkov light (DIRC) is used for particle identifica- 
tion. 

We select hadronic events using track multiplicity and 
event topology criteria. At least three reconstructed 
tracks, each with transverse momentum greater than 100 
MeV/c, are required in the laboratory polar angle region 
0.41 < #i a b < 2.54, where #i a b = is the e~-beam direc- 
tion. The event must have a total measured energy in the 
laboratory frame greater than 4.5 GeV. In order to help 
reject non-PP background, the ratio of Fox- Wolfram mo- 
ments H2/H0 is required to be less than 0.5 0. 

For reconstruction of P — > P* ± P T decays, all daugh- 
ter tracks are required to pass within 10 cm in z and 1.5 
cm in r — of the center of the beam crossing region. A 
track is identified as a charged kaon candidate using the 
Cherenkov angle measured in the DIRC and energy loss 
information (dE/dx) from the DCH and SVT. 

Neutral pion candidates are composed of pairs of pho- 
tons in the EMC. The photons must each have energy 
above 30 MeV, and their energy must sum to greater than 
200 MeV. The 7r° candidates must have an invariant mass 
between 115 and 150 MeV/c 2 . A mass-constrained fit is 
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imposed on ir° candidates, in order to improve resolution 
on the energy of reconstructed B candidates. 

We require — > tt^tt candidates to have an invari- 
ant mass within 15 MeV/c 2 of the nominal K® mass fiof . 
The transverse flight distance of the K® from the primary 
event vertex is required to be greater than 2 mm. 

To form D candidates, kaon candidates are combined 
with other tracks, assumed to be pions, and ir° can- 
didates in the event. We reconstruct D° candidates 
in the four modes K~ir + , K^tt + tt°, K^tt + tt^tt + , and 
K®ir + ir~, and D + candidates in the modes K~tt + it + 
and K®% + . We require D° and D + candidates to have 
reconstructed invariant masses within 20 MeV/c 2 of their 
respective nominal masses, except for D° decays with a 
7T° daughter, which must be within 35 MeV/c 2 of the 
nominal D° mass. Mass-constrained fits are applied to 
D° and D + candidates in order to improve the measure- 
ment of the momentum of each D. The D* + is then 
reconstructed in its decay to D°tt + . 

To select neutral B candidates from pairs of well- 
reconstructed D*^ and D T candidates, we form a likeli- 
hood that includes all measured D* and D masses: 



C = G(m DT , ov, 



H(8m D - 



core 
± ) a 8m n 



G(m D o,CTm D o) ■ 
,± ; /core): 



tail 

> "8m r 



(3) 



where the D*^- D° mass difference is denoted by 8mo*- 
Each G represents a Gaussian distribution, and H is the 
sum of two Gaussian distributions, for the core and tail 
of the 5m£,,± distribution, respectively. For a mD we use 
values individually computed for each D candidate, while 
for asm D , we use values obtained from an inclusive D* 
data sample: 0.35MeV/c 2 for the core Gaussian distri- 
bution and 1.27MeV/c 2 for the tail, and a core fraction 
(/core) of 51%. Likelihood cuts are set individually for 
each combination of D*^ and iU decay modes, using a 
detailed Monte Carlo simulation, in order to maximize 
the expected signal sensitivity. In events with more than 
one B° candidate, we choose the candidate with the high- 
est likelihood value. 

A B — > D^D^ 1 candidate is characterized by two 
kinematic variables: the beam-energy substituted mass, 

toes = *J (\fs/2) 2 — P*b 2 i an d the difference of the B can- 
didate's measured energy from the beam energy, AE = 
E* B - (y/s/2). E* B (p* B ) is the energy (momentum) of the 
B candidate in the e + e~ center-of-mass frame and ^/s is 
the total center-of-mass energy. The signal region in AE 
is defined to be \AE\ < 18MeV. According to Monte 
Carlo simulations, the width of this region corresponds 
to approximately twice the signal resolution in AE . 

The B — ► D*^D T decay candidates in the region 
5.27 < toes < 5.30GeV/c 2 and \ AE\ < 18MeV are used 
to extract CP-violating asymmetries. A sideband, de- 
fined as 5.20 < toes < 5.27GeV/c 2 and \ AE\ < 18MeV, 
and a "large sideband," defined as 5.20 < toes < 
5.27GeV/c 2 and \AE\ < 200 MeV, are used to extract 



various background parameters. The total numbers of 
selected events in the signal region, the sideband, and 
the large sideband are 197, 461, and 5187, respectively. 

To extract the number of signal events above back- 
ground, as well as the time-integrated CP asymmetry A 
(see Eq. 2), we use an unbinned extended maximum like- 
lihood fit to the Toes distribution of the £)* ± D =F candi- 
dates, including the sideband. The Toes distribution for 
the simultaneous fit to all the selected events is described 
by Gaussian distributions for the D* + D~ and D*~D + 
signals, an ARGUS threshold function 0, and a Gaus- 
sian distribution to describe a small potential "peaking" 
background contribution (due to B decays such as B° 
— > D*~D+ that are similar to the signal modes). The 
endpoint of the ARGUS function is fixed to the average 
beam energy. From studies performed with both data 
and Monte Carlo simulations, the "peaking" contribution 
is estimated to be 12 ±8 events. There are a total of four 
free parameters in the nominal fit: the shape and normal- 
ization of the background ARGUS function (2), the total 
B — > D* ± D T signal yield (1), and the CP asymmetry A 
(1). The total B — > D* ± £) T signal yield is determined to 
be 113 ± 13 events. Figure 2 shows the toes distributions 
for B -> D*-D+ and D*+ D~ candidates. 

We use a Monte Carlo simulation of the BABAR de- 
tector to determine the efficiency for reconstructing the 
B — > Z?* ± D =F signal. The efficiencies range from 6% to 
18%, depending on the D decay modes. From these effi- 
ciencies and the total number of recorded BB pairs, and 
assuming the T(45) -> B a B° branching fraction to be 
50%, we determine the branching fraction for neutral B 
to D*±D^ to be 

B(B -> D^D^) = (8.8 ± 1.0(stat.) ± 1.3(syst.)) x 10~ 4 . 
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FIG. 2: The m E s distributions of a) B -> D*~ D + and b) 
B -> D* + D~ candidates with \AE\ < 18 MeV. The fit in- 
cludes Gaussian distributions to model the signal and a small 
peaking background component, and an ARGUS function [Tl| 
to model the combinatoric background shape. 



Systematic uncertainties on the branching fraction are 
dominated by uncertainty on the charged-particle track- 
ing efficiency (8.9%), uncertainties on the branching frac- 
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tions of the D decay modes (7.4%) ^(J> an d the uncer- 
tainty on the amount of peaking background (6.8%). The 
total systematic uncertainty from all considered sources 
is 14.5%. The result is consistent with Ref. 
The fitted value for A is 

A = -0.03 ± 0.11(stat.) ± 0.05(syst.). 

Systematic uncertainties on A are dominated by poten- 
tial differences in the reconstruction efficiencies of posi- 
tively and negatively charged tracks (0.04), and by un- 
certainty in the tties resolution for B — > D* D T signal 
events (0.03). 

The method for extracting time-dependent CP asym- 
metries shares many of the techniques that are used 
for the measurement of sin2/3 in charmonium decays in 
BABAR Q . We use the same algorithms for determination 
of the flavor of the tagging B in the event, for determin- 
ing the distance Az between the B — » _D* ± Z? T and tag- 
ging B decay vertices, and for performing the maximum 
likelihood fit. We also use the same data sample, Bfl av , 
of fully reconstructed B decays to DW ± (7t =f , p T , af ) to 
measure tagging performance and Az resolution. 

The B flavor-tagging algorithm relies on the correla- 
tion between the flavor of the b quark and the particle 
types, momenta, and charges of the remaining tracks in 
the event. A multivariate algorithm is used to separate 
events into four tagging categories and to determine tag 
flavor, the details of which are given in Ref. 0. 

The elapsed time At between the B — ► D* ± D T and 
tagging B decays is determined from the measured dis- 
tance Az between the z positions of the two B decay ver- 
tices and from the known boost of the e + e~ system. A 
detailed description of the algorithm is given in Ref. 0] . 
We accept events with a^t < 2.5 ps and |Ai| < 20 ps, 
where (ta* is the error on At. We find that 93% of signal 
candidates satisfy these requirements. 

We determine the time-dependent CP asymmetry pa- 
rameters using a simultaneous unbinncd maximum likeli- 
hood fit to the At distributions of the D*^D^ and -Bfl av 
candidates, including mgs sideband samples for back- 
ground parametrization. The At distribution for D*^D^ 
signal events is described by Eq. 1. The At distribution 
of -Bfl a v events is also described by Eq. 1 with C = 1 
and 5 = where the superscript +(— ) refers to oppo- 
site (same) flavor events, comparing the reconstructed 
and tag B mesons. The mistag fraction w reduces the 
measured S and C coefficients by a factor 1 — 2u>; this 
fraction is measured within the fit for each tagging cate- 
gory, utilizing the large Bfl av sample. We convolve the At 
distribution with a resolution function modelled by the 
sum of three Gaussian distributions. The At resolution 
is dominated by the tag vertex z-position resolution and 
is parametrized in the same way as for the charmonium 
sin2/3 measurement; this is described in detail in Ref. [l3| . 
Both continuum and BB backgrounds are incorporated, 



each with a At distribution that is determined within the 
fit, using the wies sideband. 

There are 37 fitted parameters in the combined fit for 
time-dependent CP asymmetries: the CP asymmetry pa- 
rameters S_ + , C_ + , S + _, and C + _ (4); the average mistag 
fractions Wi (4), and the differences Auii between B° and 
B° mistag fractions (4), where i is one of the four tag- 
ging categories; parameters for the signal At resolution 
function (8); and parameters for background time depen- 
dence (6), At resolution (3), and mistag fractions (8). 
The -Bflav sample constrains all parameters except the 
CP asymmetries. In the nominal fit, we fix r^o = 1.542 
ps and Am d = 0.489 ps" 1 \u%. 

The time-dependent CP asymmetry fit to the B — ► 
£>*±£)=F anc j iJ flav samples yields 

S_ + = -0.24±0.69(stat.)±0.12(syst.), 

C_ + = -0.22 ± 0.37(stat.) ± 0.10(syst.), 

SV- = -0.82±0.75(stat.)±0.14(syst.), 

C + _ = -0.47±0.40(stat.)±0.12(syst.). 

The correlation between S*_ + and C_ + is 0.16 and be- 
tween S + _ and C + _ is —0.01. Besides these correlations, 
the magnitudes of all correlations of the S and C param- 
eters with any other free parameter are each less than 
0.04. Figure 3 shows the At distributions and asymme- 
tries in yields between B° and B° tags for the D*~ D + 
and D* + D~ samples, each overlaid with a projection of 
the fit result. 

Systematic uncertainties on the time-dependent CP 
asymmetry parameters are dominated by uncertainties 
in the amount, composition, and CP asymmetry of the 
background in the selected D* ± D =F events (resulting in 
errors on the parameters ranging from 0.07-0.10); the 
parametrization of the At resolution function (0.01-0.06); 
possible differences between the -Bfl av and D* ± D =F mistag 
fractions (0.01-0.04); the error on a small correction to 
the fitted asymmetries due to the limited size of the 
D* ± D =F sample (0.01-0.02); and the potential presence 
of a small amount of CP-violating interference between 
leading order and doubly-CKM-suppressed decay chan- 
nels of the tagging B meson (0.01-0.03). 

In summary, we have measured the branching fraction 
and CP- violating asymmetries for B — > D* ± I? =F decays. 
The small size of the D* ± £) =F sample currently precludes 
the observation of CP violation in this first measurement 
in this channel; however, with the addition of more data, 
future results may provide important information about 
sources of CP-violation in the B-meson system. 
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FIG. 3: Distributions of At for B -» D* + ' D~ candidates in 
the signal region with a) a B° tag (N B o ) and b) with a B° tag 
(JVgo), and c) the raw asymmetry (N B o —Ng )/(N B o + N B - ). 
The solid curves are the fit projections in At. The shaded 
regions represent the background contributions. Figures d), 
e), and f) contain the corresponding information for D*~ D + . 
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